On the occurrence of Berezinskii-Kosterlitz-Thouless behavior in highly anisotropic 

cuprate superconductors 
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The conflicting observations in the highly anisotropic Bi2Sr2CaCu2C>8+<5, evidence for BKT be- 
havior emerging from magnetization data and smeared 3D-xy behavior, stemming form the temper- 
ature dependence of the magnetic in-plane penetration depth are traced back to the rather small 
ratio, £T /f~ = ^to/^co — 0.45, between the c-axis correlation length probed above (£T) and below 
(Cc~) T c and the comparatively large anisotropy. The latter leads to critical amplitudes £jj which 
are substantially smaller than the distance between two CuC>2 double layers. In combination with 
CcVCcT — 0.45 and in contrast to the situation below T c the c-axis correlation length exceeds 
the distance between two Cu02 double layers very close to T c only. Below this narrow temperature 
regime where 3D-xy fluctuations dominate, there is then an extended temperature regime where the 
units with two Cu02 double layers are nearly uncoupled so that 2D thermal fluctuations dominate 
and BKT features are observable. 

PACS numbers: 74.72.-h, 74.25. Ha, 64.60.Fr 



Since the pioneering work of Berezinskii[lj, Kosterlitz 
and Thoulessfll (BKT) on the BKT transition in the 
two-dimensional (2D) XY model, much efforts have been 
devoted to observe the universal behavior characteris- 
tic of the KT transition, as the universal jump of the 
superfluid density @, measured in 4 He superfluid films, 
or the non-linear I - V characteristic, observed in thin 
films of conventional superconductors Signatures of 
BKT physics can be expected also in layered supercon- 
ductors with weak inter-layer coupling. Potential candi- 
dates are underdoped cuprate superconductors where the 
anisotropy increases with reduced transition temperature 
T c 0. Recent studies of the I-V characteristic 0, the fre- 
quency dependent conductivity 3, the Nernst signal Bl, 
the magnetization [lfjl [TlT| . and of the resistance [lj] have 
been interpreted as signatures of BKT behavior. On 
the other hand, several experiments ID, 14 , 15, 1(| 17, 
18, [3, 0, [5l[ failed to observe any trace of the uni- 



versal jump in the superfluid density around T c . In- 
deed, a systematic finite-size scaling analysis of in-plane 
penetration depth data taken on films and single crys- 
tals of the highly anisotropic Bi2Sr2CaCu20s+,5 reveal 
a smeared transition, consistent with an inhomogeneity 
induced finite-size effect 22]. Furthermore, there is evi- 
dence that in small thin-film samples on insulating sub- 
strates, edge effects modify the vortex- vortex interaction 
making it short-range, unlike the logarithmic long-range 
interaction needed for the BKT transition. This appears 
to make the BKT transition impossible in thin films of 
any size if they are supported by a non-superconducting 
substrate [23j. 

In this work we attempt to unravel the conflicting ob- 
servations on the highly anisotropic Bi2Sr2CaCu20s+5: 
evidence for BKT behavior emerging from the magnetiza- 
tion data of Li et al. (Io| and evidence for smeared 3D-xy 
behavior, stemming form the temperature dependence of 
the magnetic in-plane penetration depth 13,16|, 22|. Al- 



though Bi2Sr2CaCu208+a is highly anisotropic [24j, re- 
sistivity measurements [25] uncover clearly the three di- 
mensional nature of the transition in this extreme type II 
superconductor. Accordingly, sufficiently close to T c ho- 
mogeneous samples are expected to exhibit 3D-xy crit- 
ical behavior. To explore how close this should be, we 
consider the magnetic in-plane penetration data shown 
in Fig. [l] derived from the complex conductivity mea- 
surements of Osborn et al. [161 ] on epitaxially grown 
Bi2Sr2CaCu208+5 films, using a two-coil inductive tech- 
nique at zero applied. For comparison we included the 
leadin g 3 D-xy behavior given by the universal relation 



1 16TT 3 k B T 



C = Co I* 



with 



t = T/T c - 1, v ~ 2/3. 
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Apparently, there is a rounded transition pointing to a 
finite size effect, preventing the c-axis correlation length 
£~ to grow beyond the limiting length L c , set by inho- 
mogeneities. Indeed the extreme in dRe(p) /dT exhibits 
an inflection point at T = T Pc ~ 83.77 K where £~ at- 
tains the limiting length L c . As shown previously [22l ]. 
Eq. JO yields together with T c = 84 K, T Pa ~ 83.77 K, 
Re(p (T c )) /Re(p (T = 0)) = 0.011, Re(p(T = 0)) = 0.6 
K, Re(p(T p J) = 13.66 K, A afc (0) = 265 nm, and L c = 
£<T (T Pc ) the estimates 
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for the critical amplitude of the c- axis correlation length 
below T c and the limiting length L c set by inhomo- 
geneities. The rather small value of t;~ , reflecting the 
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FIG. 1: Re(p) oc 1/X 2 ab and dRe(p)/dT vs. T for the 
Bi2Sr2CaCu208+a film B, 616 A thick, derived from Osborn 
el gl.fl6l| The solid and dashed lines indicate the leading crit- 
ical behavior of a homogeneous bulk system according to Eq. 
HI where Re(p) = 650 |t| 2/3 with T c = 84 K. 



FIG. 2: J" and f+ vs. T for £~o = 2 A, £+ = 0.9 A and 
T c = 84 K, where £± = |i|~ 2/3 . The dotted line marks 
s = 15A, the distance between two CUO2 double layers. The 
arrow marks T = 83.25 K where in Fig. [T] consistency with 
3D-xy critical behavior sets in. 



0.45. 
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This ratio implies a substantial shrinkage of the 3D-xy 
fluctuation dominated regime above T c . To illustrate this 
point we included in Fig. [5] the temperature dependence 
of . Taking again > 3s = 45 A to locate the regime 
where 3D-xy fluctuations dominate we obtain roughly 
T > 83.8 K which is rather close to the T c ~ 84 K of the 
fictitious homogeneous system. 

Given the reduced 3D-xy critical regime above T c the 
system corresponds to a stack of nearly independent units 
with two Cu02 double layers as long as < 3s. In 
this intermediate regime, 2D and in particular BKT fea- 
tures are then expected. To check this expectation quan- 
titatively we consider the 3D-xy scaling expression for 
the susceptibility in the limit T > T c and H c — > 



high anisotropy, 7 = £ abo /£ c o> was confirmed in other lijl I2IL El 

Bi2Sr2CaCu208+5 films and single crystals in terms of a 
detailed finite size scaling analysis [22|. 

Although the c-axis correlation length £~ increases by 
approaching T c , the occurrence of 3D thermal fluctua- 
tions requires that in F^S^CaC^Os+a it exceeds s = 15 
A, the distance between two CuC>2 double layers. In Fig. 
[5] we depicted the temperature dependence of Even 
though it exceeds 15 A around 80 K , a glance to Fig. [1] 
shows that consistency with 3D-xy behavior is achieved 
much closer to T c only, namely above T w 83.25 K. Ac- 
cordingly, the occurrence of 3D-xy behavior requires £~ 
to exceed the distance between two CuC>2 double lay- 
ers by at least a factor of three. Above T c the situation 
is even worse because the ratio between the correlation 
lengths above (£+) and below T c is a universal quan- 
tity and in the 3D-xy universality class given by 
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where m = M/V is the magnetization per unit volume 
and Q + C^ ~ 0.9 a universal number. In type II su- 
perconductors, exposed to a magnetic field Hi in direc- 
tion i, there is also the magnetic field induced limiting 
V^o7 (aHi) with a 



3.12 




related to 
As 



length L Hi 

the average distance between vortex lines [2 
the magnetic field increases, the density of vortex lines 
becomes greater, but this cannot continue indefinitely, 
the limit is roughly set on the proximity of vortex lines 
by the overlapping of their cores. Due to these limit- 
ing length the phase transition is rounded and occurs 
smoothly. Indeed, approaching T c from above the corre- 
lation lengths combination increases but is bounded 
by L 2 Hk = $0/ (a-Hfc) where i ^ j ^ k. In this context 
it is important to recognize that the confinement effect 
of the magnetic field in direction i on fluctuations within 
a region L Bi acts only in the plane perpendicular to the 
field. Therewith Eq. © reduces for T ~ T c and H c -> 
to 
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Q + C+ k B Tj 
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because ^b/^c — 7 an d the growth of is limited by 
^ h = Lh c — y/$o/ {aH c ). On the other hand, in the 
strict 2D case f + cannot grow beyond the thickness d s 
of the system. Consequently, as —m/H c and £^~ b initially 
increase with reduced temperature in a fixed field they 
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FIG. 3: M/H vs. H for BiaSraCaC^Os+s at T = 84 K 
derived from the magnetization data of Li et al. 10] with H 
applied along the c-axis. The solid line is M(T = 84K)/H = 
6.6/ H corresponding to Eq. (0. 



then saturate due to the magnetic field induced finite size 
effect. In this case Eq. ([5) reduces to 
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In Fig. [5] we show -M/H vs. H for Bi 2 Sr 2 CaCu20 8+5 
close to T c at T = 84 K derived from the magnetization 
data of Li et al. [To] ] , with H applied along the c-axis. 
The comparison with the characteristic 2D-behavior ([7]) 
reveals remarkable agreement, extending over nearly two 
decades of the applied magnetic field. Noting that down 
to the lowest applied magnetic field, H c = 5 Oe, there 
is no sign of leveling off, arising from an inhomogene- 
ity induced finite size effect, it follows that the length 
of the homogenous regions in the afc-plane, L a t, ex- 
ceeds the attained magnetic field induced limiting length 
Lh c = \/$o/ (aH c ) = 1.14 1CT 4 cm. This behavior does 
not confirm the occurrence of intermediate 2D critical be- 
havior in sufficiently anisotropic systems above T c only. 
In addition, it uncovers that in the afo-plane excellent ho- 
mogeneity can be achieved and an applied magnetic field 
leads to the outlined finite size effect. 

To substantiate the occurrence of the magnetic field 
induced finite size effect further, we consider the temper- 
ature dependence of the magnetization at fixed magnetic 
field. When £+ < 3s and L Hc » £+ b Eq. © reduces to 
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In this regime the in-plane correlation length is ex- 
pected to exhibit the characteristic BKT-behavior in 



FIG. 4: —Mj {TH C ) vs. t = T/T c - 1 for BizSraCaCuaOg+i 
at H c = 10 Oe with T c = 84 K derived from the 
magnetization data of Li et al. QJJ]. The solid line is 
-M/(TH C ) = 10~ 8 exp(2/t 1/2 ) corresponding to Eq. © 
with f+ b oc exp(6t -1//2 ) with b = 1 and the dashed one 
-M/(TH C ) = 8.7 10" 3 . The inset shows the finite- 
size scaling function in terms of (— M/ (TH C )) / exp(26t _1//2 ) 
vs. exp(2bt~ 1/2 ) with 6 = 1. The dotted line is 
(-M/ (TH C )) / exp(2ta- 1/2 ) = 0.007/ exp(2fa- 1/2 ), indicat- 
ing the limiting magnetic field induced finite-size behavior and 
the solid one is (-M/ (TH C )) / exp(26i" 1/2 ) = 10~ 8 marking 
the BKT limit. 



the reduced temperature t — T/T c — 1, namely £+ b = 
aexp(6t _1 / 2 )[2|, with b w 1 a non- universal constant 
and a related to the vortex core diameter. However, as 
t decreases and with that £i approaches the limiting 
length Lh c = \f§o] ( a H c ) a finite size effect sets in and 
£+ & saturates in the limit L Hc << £+ b to £+ 6 = L Hc . 
To describe the resulting crossover we introduce the 

finite-size scaling function S [32| in terms of (^ab) = 

(C b ) 2 s(H c (C b ) 2 /*o) = {t ab f S ([t ab f / (aL% c 

where S ((£+) 2 / (aL%J ) - 1 for L Hc » £+, 

^(O'/KJ) - *«#hJ{&? fOT L Hc « 

and substitute £+ fc in Eq. (|8J) instead of £+ fe . In Fig. H] 
we depicted —Mj {TH C ) vs. t = T/T c - 1 at H c = 10 
Oe and the inset shows the finite-size scaling function in 
terms of (-M/(TiJ c ))/exp(26i- 1 / 2 ) vs. exp(26^ 1 / 2 ). 
It is seen that by approaching the transition temperature 
—M/ (TH C ) increases but saturates to —Mj (T C H C ) = 
8.7 10~ 3 , in agreement with the behavior shown in Fig. 
[3] Sufficiently away from T c a crossover to the BKT 
behavior, requiring L# c >> £^ and indicated by the 
solid line, can be anticipated. As the crossover extends 
over a rather extended temperature regime agreement 
with the leading BKT behavior is limited and the 3D- 
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xy critical regime is not accessible. Nevertheless, the 
finite-size scaling function reveals the flow to BKT be- 
havior is attained. To provide a consistency test we use 
M(T = 8AK)/H C = 6.6/ H c (Fig. |3J) and d s = 15 A to 
estimate a from (-M/ (TH C )) / exp(26^ 1 / 2 ) = 1CT 8 © 
with the aid of Eqs. © and ©. The result is a ~ 9.2 
A, in reasonable agreement with the estimate of Li et al. 

We have seen that the conflicting observations in the 
highly anisotropic B^Si^CaC^Os+i, evidence for BKT 
behavior emerging from the magnetization data of Li et 
al.\ld^ and smeared 3D-xy behavior, stemming form the 
temperature dependence of the magnetic in-plane pene- 
tration depth [la, [lil Eil] are a consequence of the rather 
small ratio, ^f/^T = £jo/Co — 0-45, between the c-axis 
correlation length probed above (£+) and below 
T c and the comparatively large anisotropy. The latter 
leads to critical amplitudes £ c0 which are substantially 
smaller than the distance between two CuC>2 double lay- 
ers. In combination with £j ~ 0.45 and in contrast 
to the situation below T c the c-axis correlation length 
exceeds the distance between two CuC>2 double lay- 
ers very close to T c only. Below this narrow temper- 
ature regime where 3D-xy fluctuations dominate, there 
is then an extended temperature regime left where the 
units with two CuC>2 double layers are nearly uncoupled 
so that 2D thermal fluctuations dominate and BKT fea- 
tures are observable. As this behavior relies on the uni- 
versal ratio t,t He — 0-45 and pronounced anisotropy 
it is not an artefact of Bi2Sr2CaCu20s+5, but a generic 
feature of sufficiently anisotropic cuprate superconduc- 
tors. Examples include underdoped La2- a: Sr : i;Cu04 and 
YBa2Cu307_,5 where the temperature dependence of the 
in-plane superfluid density does not reveal any trace of 
the universal jump 0, 0, 0, M, EE H EE M, E | . 
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